Casimir repulsion between metallic objects in vacuum 
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We give an example of a geometry in which two metallic objects in vacuum experience a repulsive 
Casimir force. The geometry consists of an elongated metal particle centered above a metal plate 
with a hole. We prove that this geometry has a repulsive regime using a symmetry argument and 
confirm it with numerical calculations for both perfect and realistic metals. The system does not 
support stable levitation, as the particle is unstable to displacements away from the symmetry axis. 



Introduction: The Casimir force between two parallel 
metal plates in vacuum is always attractive. A long- 
standing question is whether this is generally true for 
metallic/dielectric objects in vacuum, or whether the sign 
of the force can be changed by geometry alone. More 
precisely, can the force between non-interleaved metal- 
lic/dielectric bodies in vacuum — that is, bodies that lie 
on opposite sides of an imaginary separating plane — ever 
be repulsive? In this paper, we answer this question in 
the affirmative by showing that a small elongated metal 
particle centered above a thin metal plate with a hole, de- 
picted in Fig. [TJa) , is repelled from the plate in vacuum 
when the particle is close to the plate. The particle is 
unstable to displacements away from the symmetry axis, 
so that the system does not support stable levitation, 
consistent with the theorem of Ref. 1 , We establish our 
result using a symmetry argument for an idealized case 
and by brute-force numerical calculations for more re- 
alistic geometries and materials. We also show that this 
geometry is closely related to an unusual electrostatic sys- 
tem in which a neutral metallic object repels an electric 
dipole (in fact, one can even obtain electrostatic repul- 
sion for the case of a point charge [2 ). Anisotropic parti- 
cles are essential here; a spherical particle above a perfo- 
rated plate is always attracted, although non-monotonic 
effects in an isotropic case have been suggested for the 
null-energy condition rather than the Casimir energy [3 . 

Casimir repulsion is known to be impossible for 
Id/multilayer ^ or mirror-symmetric metal- 
lic/dielectric geometries in vacuum. Interleaved "zipper" 
geometries can combine attractive interactions to yield 
a separating "repulsive" force [6 , but the sign of the 
force is ambiguous in such geometries. (In contrast, in 
this paper the objects lie on opposite sides of a separat- 
ing z = 0+ plane and the interaction is unambiguously 
repulsive.) Repulsive forces also arise for fluid-separated 
geometries [J or magnetic [8", ^ or magnetoelectric mate- 
rials [TOl[TT]. A repulsive Casimir pressure was predicted 
within a hollow metallic sphere [121 [12 , but this is contro- 
versial as it does not correspond to a rigid-body motion, 
is intrinsically cutoff-dependent [M], and the repulsion 
disappears if the sphere is cut in half [5]. Another pro- 
posal is to use "metamaterials" formed of metals and di- 
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FIG. 1: (color online) (a) Schematic geometry achieving 
Casimir repulsion: an elongated metal particle above a thin 
metal plate with a hole. The idealized version is the limit of an 
infinitesimal particle polarizable only in the z direction, (b) 
At ^ = 0, vacuum-dipole field lines are perpendicular to the 
plate by symmetry, and so dipole fluctuations are unaffected 
by the plate (for any uj, shown here for = 0). (c) Schematic 
particle-plate interaction energy U{z) — U{oo): zero at z = 
and at z ^ oo, and attractive ior z ^ so there must be 
Casimir repulsion (negative slope) close to the plate. 

electrics arranged into complex microstructures [lOl \TE\ - 
[TT] . However, no specific metamaterial geometries that 
exhibit Casimir repulsion have been proposed, and the 
theoretical result [T indicates that repulsion in the meta- 
material limit (separations ^ microstructure) is impos- 
sible for parallel plate geometries. 

Symmetry argument: We begin by establishing repul- 
sion in an idealized geometry: an infinitesimal particle 
centered above an infinitesimally thin perfect-metal plate 
with a hole. We assume the particle is electrically polar- 
izable only in the z direction and is not magnetically 
polarizable at all (the limit of an infinitesimal metallic 
"needle") and the plate lies in the z = plane [Fig.jlja)]. 
The Casimir (-Polder) energy for such a particle at posi- 
tion X is given by [18] 

t^(x) = a.M){EA^)E.{^))iid4- (1) 

Here azz is the electric polarizability of the particle in 
the z direction and {EzEz)i^ is the mean-square z com- 
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ponent of the electric-field fluctuations at position x and 
imaginary frequency uj = This expectation value is 
evaluated in a geometry without the particle (i.e. a geom- 
etry consisting of only the perforated plate in vacuum). 
Conventionally, it is renormalized by subtracting the (for- 
mally infinite) mean-square fluctuations in vacuum. One 
way to compute the expectation value is to note that 
it is related to a classical electromagnetic Green's func- 
tion via the fluctuation- dissipation theorem. More specif- 
ically, {Ez{'x.)Ez{'x.'))^ is proportional to the electric field 
Ez{x')e~^^'^ produced by an oscillating z-directed electric 
dipole p = p^ze"*^^ at position x. 

The key idea for establishing repulsion is simple: we 
find a point x such that the classical field of an oscillating 
z- directed electric dipole at x is unaffected by the presence 
of the metallic plate with a hole. It then follows that 
C/(x) = t/(oo), implying that the energy U must vary 
non-monotonically between x and oo and hence must be 
repulsive at some intermediate points. While in most ge- 
ometries no such X exists, in the perforated plate geom- 
etry this condition is achieved by symmetry at x = 0. If 
a 2;-directed electric dipole is placed at z = in the hole, 
then the electric-field lines of the dipole in vacuum are 
already perpendicular to the plate by symmetry, as illus- 
trated in Fig. [TJb); thus, the vacuum dipole field solves 
Maxwell's equations with the correct boundary condi- 
tions in the presence of the plate, and U{z = {)) = U (oo). 
Note that this is true by symmetry at every frequency 
uj (real or imaginary), because the dipole moment p at 
2; = is antisymmetric with respect to the z = mir- 
ror plane. Intuitively, the basic point is that the electric 
dipole fluctuations of the particle do not couple to the 
plate at all when z = 0. 

For large z — that is, z much larger than the hole diam- 
eter W — the presence of the hole in the plate is negligible, 
and we must have the usual attractive Casimir-Polder in- 
teraction. So, as schematically depicted in Fig. [TJc), we 
expect the interaction energy U{z) — U (00) to be zero at 
z = 0, decrease to negative values for small 2: > (lead- 
ing to a repulsive force) then increase to zero for large 
z (leading to an attractive force). If the hole is circular, 
then by symmetry the force is purely in the z direction 
and the point of minimum U is an equilibrium position, 
stable under z perturbations; however, both the theorem 
of Ref. [1] and explicit calculations show that this equilib- 
rium point is unstable under lateral (xy) perturbations of 
the particle position. In fact, numerical calculations (not 
shown) indicate that the particle is unstable to lateral 
perturbations and tilting at all separations z. 

Electrostatics: Strictly speaking, this symmetry argu- 
ment only shows that the force must be repulsive at some 
z 0: U could conceivably have multiple oscillations. To 
definitively rule out this possibility, we rely on the ex- 
plicit numerical calculations described below. However, 
on an intuitive level, the basic behavior of the force can 
be understood from electrostatic considerations. 
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FIG. 2: (a) Schematic electrostatic interaction of a dipole 
with a neutral perforated plate (side view), depicting the 
charge density a induced on the plate. Since a is positive 
for small r and negative for large r, a can be constructed out 
of a superposition of dipoles in the z = plane, oriented radi- 
ally inward about the z axis. A simple calculation then shows 
that the interaction is repulsive for small z. (b) In contrast, 
a dipole oriented parallel to the plate is always attracted, as 
one can see from the induced charge density shown above. 

To see this, let us focus on the u = i(_ = contribution 
to the Gasimir energy ([T]); we expect the contribution 
from nonzero imaginary frequencies to be qualitatively 
similar (though this expectation can sometimes be vio- 
lated, as in the inset of Fig. |4|. The uj = contribution 
is proportional to the electrostatic energy of a z-directed 
electric dipole in the presence of a neutral metal plate. 
By the same arguments as above, such an electrostatic 
dipole must be repelled from the plate for some z > 0. 
To see this explicitly, suppose there is a static dipole at 
some position (0,0,z), and consider the induced charges 
on the plate. In the limit where the plate is infinitesi- 
mally thin, we can combine the charges on the two sides 
of the plate into a single surface charge density a. On a 
qualitative level, we expect this total charge density to 
be of the form shown in Fig. [2|a), with a positive for 
small r and negative for large r. In particular, a can be 
constructed out of a superposition of dipoles in the z = 
plane, oriented radially inward about the z axis. A sim- 
ple calculation shows that vertical force on a dipole at 
(0, 0, z) from a horizontal dipole at distance r from the z 
axis is repulsive \i r > 2z and attractive \i r <2z. Thus, 
if the hole is circular with diameter W and z < 
all the dipoles will exert a repulsive force and the total 
force is necessarily repulsive. On the other hand, when 
z ^ most of the dipoles will exert an attractive force, 
so the total force is attractive. In contrast, a dipole ori- 
ented parallel to the plate is always attracted, as one can 
see from the induced charge density schematically shown 
in Fig. [2|b). This explains why an elongated shape is 
necessary for the repulsive effect (see Fig.|4|: dipole fluc- 
tuations parallel to the plate give rise to an attractive 
Gasimir force. 

We can confirm this picture by solving this electro- 
statics problem exactly in the two-dimensional (2d) case, 
where the metal plate with a hole is replaced by a metal 
line with a gap of width W . Assuming a 2d Goulomb 
force F{r 12) = ^1^2/^12, and a z-directed dipole moment 
Pz^ we find 
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^electrostatic (^) = —Pz ' /^tt/2 i 4^2^12 * 

(2) 
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FIG. 3: (color online) Exact Gasimir force for cylinder-plate 
geometry (inset) for perfect metals (computed with BEM) 
and gold (computed with FDTD); positive (shaded) force is 
repulsive. For d < 300 nm and d > 170 nm (vertical dashed 
line), the force is unambiguously repulsive as the cylinder is 
entirely above the plate. In contrast, a perfect-metal sphere 
(diameter 60 nm) is always attracted to the plate. 



The force is indeed repulsive for small with a sign 
change occurring at z = W/2. A similar calculation for 
a ^/-directed dipole yields a uniformly attractive force. 

The repulsion in the z-directed case is quite unusual, 
even in electrostatics: in almost all cases, the electro- 
static interaction between an electric dipole and a neu- 
tral metal object is attractive, not repulsive. Indeed, on 
an intuitive level, it seems almost inevitable that a dipole 
will induce a dipole moment in the metal object oriented 
so that the force is attractive. More rigorously, one can 
prove that this interaction is attractive in several different 
limits. For example, if a dipole is very far away {z oo) 
or very close to the surface of a metal object, the inter- 
action is always attractive. One can also prove that the 
force is attractive if the metal object is replaced by a di- 
electric material with a permittivity e/eo = 1 -\- S where 
< ^ <C 1, using a perturbative expansion in S. Clearly, 
a special geometry is necessary to obtain a repulsive force 
in electrostatics, and arguably in Casimir interactions as 
well by extension to uj ^ along the imaginary frequency 
axis (see concluding remarks below). 

Numerical demonstration: Moving beyond the ideal- 
ized geometry, we expect the repulsion to be robust 
under small perturbations, such as finite particle size, 
plate thickness, and permittivity. This expectation is 
validated by the explicit calculations described below. 
We utilize two recent numerical methods, evaluating the 
Casimir force at zero temperature. First, we use a finite- 
difference time-domain (FDTD) technique that computes 
the Casimir stress tensor via the Green's function [T9( l2Q], 
with a free-software implementation [2T]. Second, we use 
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FIG. 4: (color online) Exact (BEM) 2d Gasimir force for 
ellipse-line geometry (lower inset) with perfect metals and TE 
polarization (in-plane electric field); positive force (shaded) is 
repulsive. The effects of both particle width ay and line thick- 
ness t are shown, for fixed = 0.002VK. As the ellipticity 
7 = az/tty decreases or t increases, the repulsive force di- 
minishes. Upper-left inset: frequency-resolved force F{lm(jj) 
at fixed separation d = 0.2W and t = 0^: as 7 decreases, 
attractive contributions arise from small Imuj. Upper-right 
inset: ^(Imcj) at fixed d = 0.2VF and 7 = 4: as t increases, 
attractive contributions arise from large Imcj. 

a boundary-element method (BEM) that can solve ei- 
ther for the stress tensor or directly for the Casimir en- 
ergy/force via a path-integral expression [22] . 

Figure [3] shows the Casimir force for a finite-size 
cylindrical metal particle (20 x 320 nm) above a finite- 
thickness {t = 20 nm) plate with a circular hole of di- 
ameter 1 /im, considering both perfect metals and finite- 
permittivity gold, along with the force on a perfect- 
metal sphere (diameter 60 nm) for comparison. The per- 
fect metal results were computed with BEM; the oth- 
ers were computed by our FDTD technique in cylindri- 
cal coordinates, with the gold permittivity described by 
e{uj = = 1 + ^p/^^^ where ujp = 1.37 x 10^^ rad/sec 
(the omission of the loss term, which is convenient for 
FDTD [19], does not significantly affect our results). The 
force on the sphere is always attractive, while the force 
on the cylindrical particle is repulsive for a center-center 
separation d < 300 nm. Because of the finite sizes, when 
d < 170 nm the tip of the particle intrudes into the 
hole. However, there is a range of about 130 nm for 
d > 170 nm where the force is unambiguously repulsive: 
the two objects lie on opposite sides of an imaginary sep- 
arating plane. Similar behavior is seen for an infinites- 
imally thick (t = 0) plate (Figure |3|. Somewhat sur- 
prisingly, the finite-permittivity gold exhibits a stronger 
repulsive force than the perfect-metal case of the same 
geometry; this is explained below as a consequence of 
the finite thickness of the plate. 

In order to better understand the dependence on ge- 
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ometry, we use BEM to explore the parameter space of a 
simplified 2d {yz) version of the problem: a metal ellipti- 
cal particle above a metal line with a gap of width W. We 
compute the Casimir force in this setup for perfect metals 
and 2d electromagnetism, with the standard convention 
that the electric field is in the plane, as in a "TE" mode. 
(This system is equivalent to a scalar field with Neumann 
boundary conditions on the two objects). In Figure [ij we 
explore how both the ellipticity 7 = a^/ay of the parti- 
cle and the line thickness t affect the force, for a fixed 
width W and particle length = 0.002W. As 7 ^ 1, 
the elliptical particle becomes increasing circular, and re- 
pulsion diminishes due to the attractive force associated 
with dipole fiuctuations in the y direction. We find that 
the repulsive force disappears for 7 < 1.25, when t = 0. 
Similarly, as t becomes larger, one can no longer make 
the approximation that the metal line does not affect the 
field of a z-directed dipole at d = 0, and the repulsive 
effect disappears hy t ^ O.IW for 7 = 4. 

Further insight can be gained from the contribution 
of each imaginary frequency uj = to the force for a 
fixed particle-line separation d = 0.2W (roughly maxi- 
mum repulsion). The upper-left inset to Fig.[4]shows that 
as 7 decreases, the attractive contributions first appear 
at small ^, eventually making the overall force attrac- 
tive. The right half of the inset shows that, in contrast, 
a nonzero t gives rise to attractive force contributions at 
large ^. This may explain the larger repulsive force of 
real gold compared to perfect metal in Fig. [3j the finite 
skin depth of gold cuts off the large-^ contributions, re- 
ducing the attractive effects of the finite plate thickness 
(which, in this case, dominate the attractive effects of the 
finite particle size and ellipticity). 

Concluding remarks: In this paper, we have shown 
that the sign of the Casimir force in vacuum can be 
changed by geometry alone, without "cheating" by inter- 
leaving the bodies as in Ref. 6, Consistent with Ref. 1, 
the geometry described here does not support stable levi- 
tation since the particle is unstable with respect to lateral 
{xy) translation and tilting (as shown by additional 3d 
BEM calculations). As for the question of experimen- 
tal realizations, we leave this to future work, though we 
note that one approach would be to anchor the particle 
to a substrate plate made of a low permittivity material 
using a pillar made of the same kind of material. Assum- 
ing a periodic array of such pillars and a complementary 
array of holes with a unit cell of area 10(/im)^, and esti- 
mating the force per unit cell as the maximum repulsive 
force calculated in Fig. [Sj one obtains a repulsive pres- 
sure of about 10~^ Pa. This is 2 or 3 orders of magnitude 
smaller than typical experimental sensitivities [23], but 
the repulsion could be increased further by shrinking or 
optimizing the geometry. 

This geometry was motivated by the electrostatic ana- 
logue shown in Fig. |2j where a qualitatively similar 



effect is observed. Previously, another interesting non- 
monotonic Casimir effect was also seen to have an elec- 
trostatic analogue [24]. Mathematically, the mostly non- 
oscillatory exponential decay of the Casimir-force con- 
tributions for imaginary frequencies uj = tends to 
make the total force qualitatively similar to the ^ ^ 0^ 
contribution (in fact, this similarity becomes an exact 
proportionality in the unretarded, van der Waals limit). 
This suggests that one approach for discovering "interest- 
ing" geometric Casimir effects is to first find an interest- 
ing electrostatic interaction, and then seek an analogous 
Casimir system. 
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